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Abstract 


c = chord, m 


The combined effects of blade torsion 
and dynamic inflow on the aeroelastic sta- 
bility of an elastic rotor blade in forward 
flight are studied. The governing sets of 
actuations of motion (fully nonlinear, line- 
arized, and multiblade equations) used in 
this study are derived symbolically using a 
program written in FORTRAN. Stability 
results are presented for different struc- 
tural models with and without dynamic 
inflow. The study shows that symbolic and 
numerical programs written in FORTRAN can 
be conveniently used in a complicated 
helicopter-rotor aeroelastic modeling and 
analytical process. It is observed that 
for a large number of degrees of freedom 
and for fully nonlinear models, the amount 
of data needed for the symbolic program 
increases exponentially, making it incon- 
venient to consider the multiblade equa- 
tions explicitly. However, a combination 
of symbolic and numerical programs at the 
proper stage in the derivation process 
makes the obtainment of final stability 
results an efficient and straightforward 
procedure. The symbolically generated 
equations are subsequently used to investi- 
gate the influence of elastic torsion modes 
and dynamic inflow on isolated rotor 
inplane stability in forward flight. 

Results are presented for both single- 
rotor-blade models and multiblade rotor 
systems. For both soft inplane and stiff 
inplane hingeless rotors, the elastic tor- 
sion mode significantly affects the pre- 
dicted inplane damping. Dynamic inflow 
does change the magnitude of the predicted 
damping, but the influence on damping 
trends is generally small with varying 
advance ratio or elastic coupling parameter. 


Notation 


a = lift-curve slope, 2ir/rad 

b * number of blades 


C^ o = profile drag coefficient 

c mx' c my = rotor steady pitch and roll 
moments. Eg. (11) 

^T ,C H' C W = rotor steady thrust, drag force, 
and weight coefficient, Eq. (11) 

C T ,C M .G L = harmonic perturbation coeffi- 
cients of thrust, pitching 
moment, and rolling moment, 

Eq. (3) 

D = partial derivative matrix, 

Eq. (14) 

f = flat-plate area 

F = forcing function, Eq. (6) 

J = number of points used in har- 

monic analysis, Eq. (5) 

K a - blade cross-section polar radius 

of gyration, m 

K m * blade cross-section mass radius 

of gyration, m 


Kjru'Km 2 = Principal mass radii of gyration, 

m 


L = number of harmonics used in the 

harmonic analysis, Eq. (10) 

[m],[£] = dynamic inflow matrices 

[M] , [C] , ■ constant mass, damping, and 
[K] stiffness matrices, Eq. (6) 

n = number of the harmonics in har- 

monic analysis, Eq. (9) 

N * total number of blade modes used 

q * perturbation degrees of freedom, 

Eq. (15) 
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^o'Sc'^s a vectors of collective and cyclic 
modes , respectively 

R = rotor radius, m 
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Introduction 


u = dynamic inflow quantities, 

Eq. (19) 

U T' U P = tan 9 ent ^ al and perpendicular 

velocity components, m/sec 

v ^ = induced velocity 

w,v,$ = flap, lag, and torsion 

deflections 

w,v,$ = steady flap, lag, and torsion 

deflections, Eq. (2) 

x * blade coordinate along the 

radius 

z * first-order variable degrees of 

freedom, Eq. (15) 

a = rotor-shaft plane angle of 

attack, Eq. (11) 

a R * wake skew angle, Eq. (4) 

e pc = precone angle, rad 

Y = Lock number 

Aw,Av,A <J> *= perturbation flap, lag, and tor- 

sion deflections, Eq. (2) 

^K' W K * real and imaginary parts of the 
characteristic exponent 

n i'^i' 0 i “ mcx * e shapes for flap, lag, and 
torsion 

8 ■ pitch angle, rad 

Aj< ■ characteristic exponent 

X = steady inflow (free stream plus 

induced flow) 

U ■ advance ratio 

v ■ inflow parameter, Eq. (5) 

v o' v c' v s * uniform, longitudinal, and lat- 

eral inflow components 

a « solidity ratio ■ bc/irR 

ip ■ azimuth angle, nondimensional 

time 

u w *u v tO)± - nondimensional rotating flap, 
lag, and torsional frequencies 

ft - blade rotational speed, rad/sec 

(") ■ nondimeneionalized quantity, 

equilibrium deflection 

(*) ■ time derivative 


Hingeless rotor blades are less com- 
plex mechanically and provide more rotor 
control power and damping than articulated 
rotor blades. However, the complex aero- 
elastic behavior of hingeless rotors 
requires a rigorous analysis for an effec- 
tive design procedure. The modeling 
requirements of hingeless-rotor-blade 
aeroelasticity have been studied for many 
years and are briefly reviewed here. 

Initial analyses focused on the inves- 
tigation of flap-lag stability of torsion- 
ally rigid blades with spring-restrained 
hinges at the hub to simulate bending 
flexibility. The stability of this type 
of model was analyzed for both hover 1 and 
forward flight. 2 Flap-lag stability of 
elastic blades with uniform properties was 
studied by Ormiston and Hodges, 1 based on 
a derivation of nonlinear partial differ- 
ential equations suitable for elastic 
hingeless blades. Similar equations were 
studied by Friedmann and Tong. 3 Efforts 
were also made to investigate the complete 
blade problem by including blade torsional 
deflections. Friedmann and Tong 3 approxi- 
mated the torsional deflection by rigid- 
body pitching motion (root torsion) ; they 
found that torsion motion was important 
and that the stability characteristics 
were sensitive to the number and type of 
assumed bending-mode shapes used. Flap- 
lag structural coupling was not included. 
Hodges and Ormiston 11 presented extensive 
numerical results for the stability char- 
acteristics of elastic hingeless blades 
with flap-lag-torsion motion in hover. 

They found that torsional deflections of 
hingeless rotor blades are strongly influ- 
enced by the nonlinear structural moments 
caused by flap and lead-lag bending. This 
bending-torsion structural coupling is 
proportional to the product of the flap 
and lead-lag bending curvatures and to the 
difference between the two bending flexi- 
bilities. This study also showed the 
effect of precone, structural coupling, 
and torsional rigidity on the isolated 
blade stability boundaries. 

Friedmann and Kottapalli 5 analyzed the 
coupled flap-lag-torsional dynamics of 
hingeless rotor blades in forward flight. 
They noted that nonlinearities are impor- 
tant in an aeroelastic stability analysis 
and that forward flight is strongly coupled 
with the trim stAte. However, only flap- 
ping motion was used in calculating the 
rotor trim condition. It was observed 
that forward flight (increasing advance 
ratio) is stabilising for soft inplane 
rotors and destabilizing for stiff inplane 
rotors. In all these studies, the aero- 
dynamic forces were obtained from strip 
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j theory based on a quasi-static approxima- erning equations of motion for an elastic 

‘} tion of two-dimensional , unsteady airfoil rotor blade in forward flight. A Lagran- 

; theory. gian formulation is used to obtain the 

j equations in generalized coordinates. The 

Simultaneous efforts have been made to program generates the steady-state and 
improve the aerodynamic model used in these linearized perturbation equations in sym- 

analyses by including unsteady airflow bolic form and then codes them into 

effects. One approach is to model the FORTRAN subroutines. Subsequently the 

induced velocity as a time -dependent , coefficients for each equation and for 

| three-degree-of- freedom system. This each mode are identified through a numcri- 

: dynamic inflow theory has been applied to cal program. The harmonic balance equa- 

! rigid-blade flap-lag analyses, both in tions , if required in the calculation of 

hover and forward flight, and to the the deflected equilibrium position of the 

coupled rotor-fuselage problem in hover. 9 ' 10 blade, can also be obtained from the sym- 
It was observed that the dynamic inflow bolic program. The governing multiblade 

increased the lag regres sing-mode damping equations are derived explicitly using 
and reduced the body pitch and roll damping HESL. This is the first Lime that multi- 
for the parameters considered. These ana- blade equations are derived explicitly 

lytical results correlated well with exper- using this symbolic formulation approach 

imental results. 10 However, the conclu- to study the stability of an elastic rotor 

sions presented in Refs. 6-10 were based blade in forward flight. The multiblade 

on several restrictive assumptions; for equations are capable of accommodating any 

example, zero elastic coupling, fixed number of elastic blade modes. Because 

solidity ratio, and rigid flap-lag rotor- the complete analytical process, from 

blade models with no torsional flexibility. derivation to numerical calculation, is 

The effects of dynamic inflow and torsion automated, it is an efficient and accurate 

flexibility on the aeroelastic stability means for analyzing helicopter rotor 

of an elastic rotor blade in hover to a aeroelasticity. 

number of parameters was recently pre- 
sented. 11 It was shown that for torsion- The present study differs from previ- 

ally flexible blades, the dynamic inflow ous ones in the following respects: 

effects depend on the elastic coupling 1) symbolic manipulation with FORTRAN is 

parameter. For certain values of elastic used to derive the governing equations in 

coupling, the dynamic inflow effect may in forward flight for an elastic rotor blade? 

fact be negligible. 2) complete elastic flap-lag-torsion blade 

degrees of freedom are used for the trim 
In summary, general nonlinear differ- calculation? 3) explicit multiblade equa- 

ential equations for the elastic rotor tions are derived symbolically for stabil- 

blades (used in the above analyses) have ity calculation to compare with the single- 

been developed by several researchers. 12 “ 14 blade solution? 4) dynamic inflow is 

These models have elastic flap, lead-lag, included in the aeroelastic stability 

and torsion degrees of freedom, with non- solution of an elastic blade in forward 

linearities owing to moderate elastic flight? and 5) damping data in forward 

deflections, in those studies, it was flight are presented for varying elastic 

observed that for a given ordering scheme, structural coupling, 

the final equations differed by a number of 

small nonlinear terms. These differences To demonstrate the usefulness of this 

depend in part on the stage at which, in analytical capability, stability results 

the process of derivation, the ordering are presented for several hingeless-rotor- 

scheme is applied. The application of the blade structural models. The influence of 

ordering scheme at a later stage in the dynamic inflow in forward flight with an 

derivation process requires much time in elastic hingeless rotor is also investi- 

deriving and independently checking the gated. The hingeless-rotor stability 

final equations. This has led to attempts results presented in this paper using the 

to share the algebra with computers through symbolic program reflect the combined 
symbolic processors. Both general and effect of an improved structural model 

special purpose programs have been devel- (by including torsion) and an improved 

oped and are available. 18 The program aerodynamic model (by including dynamic 

} ' i copter Equations for Stability and inflow) . Results are presented for elas- 

Locis (HESL) , appropriate to rotary-wing tic blade flap-lag-torsion analysis and 

aeroelasticity investigations written in for flap-lag analysis with and without 

FORTRAN IV, was presented in Ref. 16. The dynamic inflow, 
approach used in developing this program 
and its use in analyzing the aeroelastic 

stability of an elastic rotor blade in Formulation 

hover was presented in Ref. 11. 

„ _ Figure 1 shows an elastic blade with 

In the present paper, the HESL pro- the coordinate system used in this study, 

gram has been extended to derive the gov- The blade has uniform mass and stiffness, 
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no twist, and no chordwise offsets of the 
elastic axis, tension axis, or center of 
mass. The elastic axis is coincident with 
the x-axis of the x,y,z coordinate sys- 
tem rotating with a constant angular veloc- 
ity (ft) about a fixed point at the origin. 
The y-axis lies in the plane of rotation, 
and the x-axis is rotated through a small 
angle (Bp C ) from the plane of rotation. 

The deflections of the beam are u (axial 
deflection) , v (lagwise bending) , and w 
(flapwise bending) of the elastic axis 
parallel to the x,y,z coordinates, 
respectively. A second coordinate system, 
x 1 , y 1 , and z', is fixed to the blade, with 
y f and z' axes parallel to the principal 
axes of the deformed blade cross section. 
This coordinate system moves with the blade 
cross section as it undergoes bending, tor- 
sional displacements, and pitch angle (6) 
rotation. Before deformation, the princi- 
pal axes of the blade cross section are 
rotated with respect to the undeformed 
coordinates by the pitch angle. After 
deformation, the elastic axis is displaced 
by u,v,w, and the blade is twisted through 
the angle <f>. The aerodynamic inflow 
dynamics couple with the blade dynamics as 
a feedback loop (Fig. 2) . The total inflow 
(vj_) is assumed to consist of a steady 
value (X) and dynamic inflow components 
(v 0 , v c , and v 8 ) that vary with time. 

In this study, the entire problem 
formulation is performed by the computer; 
there is minimum user interface other than 
specifying blade geometry and the desired 
blade model representation. In general, 
the formulation of the rotary-wing aero- 
elastic problem consists of the following: 
writing the transformation matrices between 
the coordinate systems before and after 
deformation; calculating the position vec- 
tor of a mass point of the deformed blade 
section; forming strain displacement rela- 
tions; and calculating stresses and air 
velocity components in the flap, lag, and 
torsion directions (see Refs. 12-14 for 
more details) . These expressions include 
geometrical nonlinearities owing to the 
assumption of small strains and moderate 
deflections which give rise to numerous 
higher-order nonlinear terms. So an order- 
ing scheme, based on assigning orders of 
magnitude to the various physical param- 
eters, is used to reduce the number of 
terms. The governing equations of motion 
are then obtained using Hamilton's princi- 
ple. These equations are nonlinear, par- 
tial differential equations in u,v,w, and 
0 deflections. These are converted to 
ordinary differential equations using 
Galerkin's method by expressing the bend- 
ing and torsion deflection* in terms of 
generalised coordinates and mode-shape 
functions , 


n f 

W = 53 Rw i (<P) Hi (x) 

i=i 

N L 

V - E RViOMCiW (la) 

i=i 

Ny 

$ = W &i (x) 

i=i 

and by expressing the induced velocity as 
v^ * X + v Q + v c R cos if; + v s x sin i|» 

(lb) 

where ip * ftt, S * x/r, and ni,Si,0i are 
mode shapes; R is the blade radius; and 
N F , N l , and Nip are the numbers of flap, 
lag, and torsion modes, respectively, used 
in the analysis. In this study uncoupled 
rotating modes evaluated st zero pitch are 
used. This yields N nonlinear, non- 
homogeneous ordinary differential equa- 
tions in terms of modal generalized coor- 
dinates Wi, vj,, and 4>i , where N is the 
total number of flap, lag, and torsion 
modes used in the analysis. The equations 
have periodic coefficients in the mass, 
damping, and stiffness matrices. These 
equations are then linearized for small 
perturbation motions about the deformed 
blade time-dependent equilibrium position 
by expressing the generalized coordinates 
in terms of the equilibrium quantities and 
small perturbation quantities: 

w i “ ^i(v) + Awi(ip) 

Vj^ = (ip) + A Vi (if;) (2) 

4>i * ii(^) + A$i (ifi) 

Two sets of equations are obtained 
from this operation? a set of N_ non- 
linear equations in wi, v*, and $i, which 
define the deflected equilibrium position 
of the blade, and a set of N equations 
obtained by subtracting the equilibrium 
equations and discarding all nonlinear 
products of the perturbation quantities, 
AWi, Avi, and A$i. Three more equations 
are obtained for the dynamic inflow com- 
ponents from rotor perturbations in aero- 
dynamic thrust (Crp) and in pitch (C^) and 
roll (Cx,) moments (see Dynamic Inflow, 
below) , The coefficients of these equa- 
tions are also functions of the equilib- 
rium solution. 
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Dynamic Inflow 

The dynamic inflow equations are 
related to the blade degrees of freedom 
through the variations in thrust, pitch- 
ing, and rolling moments: 



P°1 

H°1 


r cT i 

[m] < 

v c V + 


► = 4 



l v sJ 

UsJ 


UJ 


The elements of [m] and (£] define the 
various dynamic inflow models that can be 
included in an analysis. Reference 8 pre- 
sents a hierarchy of models having differ- 
ent elements for (m] and [£] from actuator 
disk theory in forward flight. The ele- 
ments of [£] depend on the wake skew angle 
at the rotor: 

°R = tan" 1 [A] (4) 

where A is the steady inflow. 

Of the 13 models presented in Ref. 8, 
the partially constrained model gave good 
results. In the present paper, this par- 
tially constrained theory is used to 
obtain the dynamic inflow results. The 
elements of [m] and [£] are given by 


^28 
11 " 75ir 


nij j ® m. 


-16 


2 2 ~ ‘"3 3 - 457, ' 

n»ij = 0 , i # j 


£ = I » .15* ( j ~ sin a R V /2 

n 2 ' is 6T \1 + sin a ' 


'22 1 + sin a B 


£ ** £ 

31 13 


-4 sm a* 


53 1 + sin ' 

*12 * *21 = *23 * *32 = 0 


£j j - It.., v - t an oj 

ij v xd ' + x*)W* 


(5) 


Equations from HESL 

The governing equations of motion of 
the rotor blade are derived using HESL with 
two modes for each blade degree of freedom. 
The principles involved in the development 
of the symbolic program HESL are described 
in detail in Refs. 11, 16, and 17. Basi- 


cally, the program assigns numbers to the 
variables forming the required expressions 
and then manipulates those numbers to 
obtain the required algebraic quantities. 
The integration, differentiation, pertur- 
bation, and multiblade coordinate trans- 
formation are performed by substituting 
known relations required for these opera- 
tions. The symbolic program can handle 
both individual expressions and matrices. 
The program generates the steady and per- 
turbed equations in a single operation and 
outputs them individually. This is con- 
venient in the case of forward flight 
because of the large number of terms pres- 
ent in each equation; it is also convenient 
for the different analytical processes 
required for the steady and perturbed sets 
of equations. 

The inputs to the program are the 
relations, in alphanumeric format, for the 
position vector, for the strain expres- 
sions, for the air-velocity components, 
and for the transformation matrices as 
given by Kaza and Kvaternik. 13 The inte- 
gration relations (if known) , differentia- 
tion relations, the order of the variables, 
the ordering scheme to be used, and the 
variables for which coefficients are to be 
collected are also given as data. In the 
present paper, the order of the variables 
and the ordering scheme used are the same 
as those followed in Ref. 4. All the 
0(e) 2 terms, compared to 0(1), except 
those that contribute to lead-lag and tor- 
sion damping, are neglected. Nonlinear 
rate products (vw, v 2 , etc.) are retained 
since they contribute to the linearized 
stability analysis. Although any general 
ordering scheme could have been used to 
obtain the final equations of motion, this 
ordering scheme is considered representa- 
tive and adequate for demonstrating the 
capability of the symbolic analysis pro- 
cess. The program calculates the strain 
energy, kinetic energy, and generalized 
forces for a given orderin; scheme in gen- 
eralized coordinates using Eq. (1) . The 
perturbation relations as given in Eq. (2) 
are substituted to obtain the steady and 
perturbed terms. The program generates 
both the steady-state (nonlinear) and 
linearized perturbation equations and the 
loading terms necessary for an aeroelastic 
stability and response analysis. The 
rotor-thrust, pitch-moment, and roll- 
moment equations required in the dynamic 
inflow equations are also obtained using 
the perturbed aerodynamic forces. The 
equations are written into FORTRAN sub- 
routines for subsequent numerical calcula- 
tions. A numerical program subsequently 
identifies the mass, damping, stiffness, 
and forcing coefficients for each general- 
ized degree of freedom. For the results 
presented here, it took about 300 sec to 
symbolically derive both the structural 
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and aerodynamic equations on a VAX 11/780 
computer, 

A brief description of the program 
input and output follows. Figure 3 is a 
flow diagram of the aeroelastic arfalysis 
using the symbolic and numerical programs. 
Table 1 shows the FORTRAN symbol defini- 
tions used for the original variables. 

Table 2 shows the input required to calcu- 
late tangential and perpendicular blade 
cross-section velocities U T and Up/ using 
the transformation matrix LAFP (READ MATRIX) 
and the air-velocity vector VEL(READ 
MATRIX) . By multiplying the two matrices 
{FORM MATRIX) with ordering scheme *E2Dl, 
the vector: AVEL is obtained, which gives 
the components of the velocities in radial, 
tangential, and perpendicular directions. 
The vector components are redefined as 
expressions by command MATRIX EXPRESSION. 
The actual velocity components are the neg- 
ative of the original expression, and are 
therefore negated by calling the NEGATE 
command, thus giving the actual velocity 
expressions. This procedure is slightly 
different from the one presented in 
Ref. 11, where manipulations were per- 
formed at the expression level. Here the 
manipulations are extended to include 
matrix operations. It should be noted 
that for a hingeless rotor, the axial dis- 
placement can be solved for a priori as a 
function of flap and lag bending. In the 
present paper, expressions for axial dis- 
placement and axial velocity are taken 
from Ref. 4 and supplied as data to the 
program. 


Trim and Periodic Equilibrium Solution 

The nonlinear periodic coefficient 
equations obtained earlier can be solved 
for the periodic response in the time 
domain using a Floquet method or in the 
frequency domain using a harmonic balance 
method. Either will yield the time- 
dependent equilibrium position about which 
the nonlinear equations can be linearized 
for an eigensolution. In forward flight, 
this equilibrium position is coupled with 
the entire trim state of the helicopter. 

The trim state is the steady-state condi- 
tion achieved by the system as time 
increases without bound, with the controls 
fixed and no external output. Calculation 
of trim position requires establishing the 
control settings for a given flight condi- 
tion. The control settings are collective 
pitch, longitudinal and lateral cyclic 
pitch, and the rotor-shaft angle of attack. 
The induced velocity, which depends on the 
generated thrust and advance ratio, is also 
calculated. 

In this paper, the harmonic analysis 
method, coupled with an iteration on the 
trim state, is used to calculate the equi- 


librium position and the trim settings. 
This procedure 18 consists of the follow- 
ing. An iterative inner loop in which the 
solution for the periodic motion is 
obtained with fixed values of the trim 
variables is nested within an outer loop 
in which the solution for the trim vari- 
ables is obtained. The rotorcraft motion 
is solved for the periodic motion by the 
harmonic analysis method, which directly 
calculates the harmonics of a Fourier 
series representation of the motion. The 
procedure advances the rotor around the 
azimuth, calculating the forcing functions 
in the time domain and then updating the 
harmonica of the response. The forces and 
moments acting on the rotor are calculated 
from this response and the controls are 
adjusted until the equilibrium of forces 
and moments required for the specified 
operating state is achieved. 

For the harmonic analysis method, the 
governing equations of motion are written, 
with all the time-dependent and nonlinear 
terms as a forcing function, as 

[M] { X } + [C]{X> + [K] {X } = F (X,X, ip) (6) 

where matrices M, C, and K are the con- 
stant mass, damping, and stiffness 
matrices and X is the vector of degrees 
of freedom. The function F is evaluated 
at J points around the rotor azimuth 

Fj * FOPj) (7) 

and the harmonics of a complex Fourier 
series representation of F are given by 

?n ■ J £ ^ e‘ in *j (8) 

j*i 

Then the nth harmonic of the motion is 
given by 

x n * H n lp n ( 9 ) 

where 

H * K - n 2 M + iCn 

The iterative solution proceeds as follows. 
At a given azimuth i|m , the blade motion 
is calculated using the current estimates 
of the harmonics: 

X <■ £ X n eini(,j 

n— L 

( 10 ) 

X - £ X n in e ini(/ j 

n— L 

where L Is the number of terms used In 
the complex expansion of X. The forcing 
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function is evaluated with this motion. 
Then the change in the harmonics owing to 
the change in the forcing function is cal- 
culated and added to the harmonics calcu- 
lated in the last revolution. After every 
revolution, the old and new harmonics are 
checked until convergence. 


D = 3M . I - . . . J* 

U 3v [ 3 Vi 




M (v^ ) - - Av i ) 

AVj_ 



After obtaining the harmonics of the 
blade motion, the equilibrium uf the 
forces and moments is checked. If equi- 
librium is not satisfied, the trim set- 
tings are increased and the procedure is 
repeated until equilibrium is met. The 
following assumptions are made in arriving 
at the equilibrium of the forces and 
moments. The helicopter is in straight 
and level steady flight; the rotor-hub 
moments are trimmed to zero; and tail, 
fuselage moments, and side-force compo- 
nents are neglected. Then the equilibrium 
forces and moments are given by 

C T cos a + C H sin a = C w 
-C T sin 1 + C H cos a = -C x (11) 

c mx = C my = 0.0 

where C x = (1/2) fy 2 and a is the angle 
of attack of the shaft. 

In calculating thrust, horizontal 
force , and the hub moments , the steady 
inflow appears as a parameter that in turn 
depends on the rotor thrust and shaft 
angle of attack. In this paper, the 
steady inflow is taken as an equation of 
constraint and solved along with the four 
equations above ; 

X - U tan a - C t /[ 2(X* + u 2 ) 1 ^ 2 ] = 0.0 

( 12 > 

The increments in the trim settings are 
calculated using a modified Newton-Raphson 
method. 1 * If v is the control variable 
and M is the target to be achieved, then 
a first-order approximation of M(v) is 

M (target) = M n+1 - M n + jg (v„ +l - v„) (13) 


where v^ is the ith control variable 
and Avi is its increment. The conver- 
gence is checked when the tolerance level, 
as specified by the parameter e, is met. 
For more details on this procedure see 
Ref. 18. In the present paper, all the 
degrees of freedom - blade flap, lead-lag, 
and torsion degrees of freedom -are used 
in calculating the trim state and periodic 
response. 

The periodic response can also be 
obtained using Floquet methods. 3 A review 
of the use of these methods to obtain the 
initial conditions, forced response, and 
stability data is given in Ref. 19. 


Aeroelastic Stability Solution: 
Single Blade 

Once the time-dependent equilibrium 
position is determined, the nonlinear 
equations are perturbed about this equi- 
librium position, as given by Eq. (2). As 
mentioned earlier, the symbolic program 
generates the perturbed equations and 
writes them into subroutines. It should 
be noted that squares of the perturbation 
quantities are neglected by the symbolic 
program by employing the ordering scheme. 
The final system of equations for stabil- 
ity is 

[P] (Aq) + [Q] {Aq} + (R]Uq) = 0 (15) 
or 

{2} = {Ag]{z} 

The stability of this linearized sys- 
tem is determined from Floquet theory by 
evaluating the characteristic exponents 
of [A s ], 

*K s + (16) 


or 

v n +i ■ v n + D” 1 [M (target) - M n )F 

where v n and v n+1 indicate the current 
and new estimates of v, respectively, and 
F < 1 is included to avoid overshoot 
oscillations in the trim iteration by 
reducing the step size. The partial deriv- 
ative matrix D is 


The mass, damping, and stiffness terms are 
identified by a numerical program for each 
time-step. The linearized system is stable 
when Ck < 0. 


Multiblade Coordinate Transformation 

To provide a better understanding and 
to include dynamic inflow effects, which 
are referenced to the fixed system, it is 
necessary to convert the equations into a 
fixed coordinate system. In the case of 
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hovering flight, this transformation can be 
easily performed by rearranging the coeffi- 
cients of the equations, since the coeffi- 
cients are constant. For forward flight, 
the degrees of freedom, as well as the 
coefficients, are periodic. Hence, the 
multiblade coordinate transformation (MCT) 
is more complicated. The MCT or Fourier 
coordinate transformation is a linear 
transformation of the degrees of freedom 
from the rotating to nonrotating frame. 

Let X be a degree of freedom {dimension- 
less) in the rot . ting frame for the ith 
blade. Then, for a three-bladed rotor, 
the relations 

= X Q + X c cos + X s sin ^ 

Xi » x 0 + (x c + Xg) cos i + (Xg - X c ) sin ^ 

* X 0 + (X c + 2Xg - X c ) cos 

+ (X g - 2X C - X c )sin ^ (17) 

give the ith blade degree of freedom, 
using multiblade coordinates in the non- 
rotating frame. The variables Xg, X c , 
and X s are the rotor degrees of freedom, 
and describe the motion of the rotor as a 
whole in the nonrotating frame. 

The MCT involves the following 
steps 20 : 1) expansion of each degree of 

freedom into multiblade coordinates ; 

2) multiplying the resulting expression 
with multiblade functions like 1, cos ip, 
sin ip, cos 2ip, sin 2ip, etc., depending on 
the number of blades; 3) replacing products 
of sines and cosines as sums of sines and 
cosines, using trigonometric relations; and 
4) deleting terms that are not multiples of 
the number of blades (summation rules) . 
Using the symbolic program, this is 
achieved as follows. The multiblade expan- 
sions of each degree of freedom and their 
time derivatives are given as a table of 
relations (Eq. (17)). The multiblade func- 
tions like 1, cos ip, sin ip, cos 2ip, sin ?ip, 
etc., require 1 in transforming the equa- 
tions, are read as data (or can be gener- 
ated within the program) . The trigonomet- 
ric relations givinq the product of sines 
and cosines as sums of sines and cosines 
are given as a table of relations. These 
require that the equation derivation be 
several runs until all the required rela- 
tions are included in the table of rela- 
tions. Then the command PERFORM MULTIBLADE 
TRANSFORMATION multiplies the equation with 
each of the multiblade functions, substi- 
tutes the multiblade expansion for each 
degree of freedom, substitutes the trigo- 
nometric relations (from the given tables 
of relations), and checks for the multiples 
of the blade harmonics. Only terms con- 
taining multiples of the number of blades 
are retained. The interface routines sep- 
arate the terms into constant and periodic 


parttf, with the coefficient of each har- 
monic separated, and writes them into sub- 
routines. In this manner, the constant- 
coefficient approximation is easily done. 
In the present paper, the perturbed gov- 
erning equations of motion and perturba- 
tions in thrust and moment equations are 
converted into multiblade coordinates. 

The multiblade solution was checked for 
a cc racy with a single-blade solution 
without dynamic inflow. It should be 
noted that the trim-value harmonics enter- 
ing as nonlinear contributions should be 
defined as symbolic data. The input data 
increase as more nonlinear terms are taken 
into account in addition to the data given 
for the multiblade expansion of the 
degrees of freedom. However, the output 
may be smaller since only terms that are 
multiples of the number of blades are 
retained. Since this is a feasibility 
study undertaken to obtain explicit multi- 
blade equations using a symbolic program 
in FORTRAN, the nonlinear quantities are 
assumed to provide only first harmonic 
forcing contributions. For the results 
presented here, the program was rur. on a 
VAX 11/780 computer. It took about 
250 sec to derive tb'j multiblade equations 
for each blade decree of freedom, and 
about 120 sec to write these into subrou- 
tines for numerical analysis. 

It should be noted that by giving the 
expansion of each degree of freedom into 
its harmonics and by giving the trigono- 
metric relations as data to the symbolic 
program, explicit harmonic balance equa- 
tions can also be derived. However, 
because of the amount of input required 
to perform a symbolic formulation of the 
harmonic balance and multiblade equations, 
the program HESL is convenient for 
explicitly considering the symbolically 
derived equations only if the number of 
degrees of freedom is small. As pointed 
out in Ref. 20, numerical schemes are 
bettor suited to yeneral models for effi- 
ciently obtaining the harmonic balance 
equations and nultibiade equatiors after 
the steady and perturbed equations are 
obtained from the symbolic program. 

The symbolic program separates the 
terms containing the periodic variable 
cos Kc and sin Nt and writes the equation 

as 

A ( t) - A 0 + A n cos Nt + B n sin Nt (18) 

In subroutine form they are referred to a» 
A (1) , A (2) , and A(3) • This allows for 
direct elimination of the matrices A(2) 
and A (3) for a simple constant-coefficient 
approximation analysis. 
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Stability Solution: Multiblade anc 

Dynamic Inflow 

The final governing equations of 
motion can be written as 

[P] {q} + CQHq} + [RHq) + [T] {u} = 0 

(19a) 

for the blade equations, and as 
[A] {q} + [B ) {q} + [C]{q} + [G]{u} 

* [m]{il} + ur l {u} (19b) 
for the dynamic inflow equations, where 
{q} is (q c ,q c ,q g } 

{ q Q } is vector of all collective modes 

{ q c } is vector of all lateral cyclic modes 

{q g } is vector of all longitudinal cyclic 
modes 

{u} is {v 0 ,v c ,v g } 


Defining A - {u}, Eqs. (19a) and (19b) can 
be combined as 



(G 


l )j 


|{X> 


+ 


R 

C 


fx} * 0 


where {X} is 



( 20 ) 


The final stability equations in state 
vector form are 


(Y) » [A g ] { Y} (21) 

where 


(Y) 



The stability results are obtained by cal- 
culating the eigenvalues of [a 8 ] , in a 
manner similar to that used in the single- 
blade case. 


The size of the state matrix depends 
on the number of modes and blades. For the 
flap-lag-torsion model with two modes each, 
the size of the state matrix is 36 * 36 
without dynamic inf ow and 39 * 39 with 
dynamic inflow. The corresponding values 


for a flap-lag model are 24 x 24 and 
27 x 27, respectively. 


Results and Discussion 


Results are presented for a uniform 
blade with zero built-in twist, zero pre- 
ccne angle and zero blade offsets. 
Reversed-flow effects are neglected. A 
three-bladed rotor is considered. Two 
rotating modes for each flap, lag, and 
torsion degrees of freedom are used in the 
calculation. These modes are calculated 
at zero pitch and are obtained from five 
nonrotating modes. Results for both a 
single-blade solution and a multiblade 
solution are presented for different blade 
structural models. All results are for 
a propulsive trim condition, specified 
for a weight coefficient of C w /o = 0.07 
and an equivalent drag area 
D/q - f = 0 . r ~(*R 2 ), where D is the 
drag force, . 1 q is the dynamic 

pressure. 

In the derivation of the equations, 
the order of magnitude assigned for each 
parameter is the same as that followed in 
Ref. 4. The other parameters used for the 
numerical study are 

<D V = 0.7, 1.4; = 1.15; = 3.0; 

c/R = 0.07854; y = 5.*.. 0 * 0.1; 

a - 2tt; C do = 0.01; 0 pc - 0.0; 

Knu/Km* “ O* 0 '* V* * °- 025 > 

( K A /K m )* * 1.5 

I,ead-3ag damping values (real part of 
the characteristic exponent) are presented 
for a soft inplane and a stiff inplane 
rotor with and without dynamic inflow. 

The results are presented for investigat- 
ing 1) the effect of degrees of freedom 
used in the f-rim analysis on the lead-lag 
damping, 2) the effect of using only one 
torsion mode, 3) the inclusion of a 
dynamic inflow model, and 4) the differ- 
ence between periodic and a constant- 
coefficient approximation. 

Single-Blade Results 

The effect of the number of degrees 
of freedom used in the trim analysis on 
the lead-lag damping is shown in Figs. 4 
and 5. Figure 4 shows the lead-lag damp- 
ing plotted versus advance rati • for a 
soft inplane rotor (w v - 0.7). It can be 
seen that a flap-lag-torsion stability 
analysis from a flap-trim analysis under- 
predicts the lead-lag damping. The second 
mode shows the same trend with the differ- 
ence in predicted damping increasing with 
advance ratio. Figure 5 shows the lead- 
lag damping plotted for a stiff inplane 
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rotor (wy - 1.4) as a function of advance 
ratio. The results also show an increase 
in damping when a fJ ap-lag-torsion trim 
analysis is used. It Is also noted that 
for an advance ratio of 0.37 <p< 0.41, 
the first-mode roots separate and one root 
becomes less stable and the other becomes 
more stable. The damping does reduce as 
the advance ratio is increased. The second 
mode remains stable at all advance ratios 
considered. 

The increase in damping observed above 
for both soft inplane and stiff inplane 
rotors is perhaps a result of the differ- 
ent time-dependent equilibrium positions 
used. A full flap-lag-torsion trim analy- 
sis is consistent in that the blade model 
has the same degree of complexity in both 
the trim ana the stability analysis. It 
should be noted that qualitatively the same 
type of trend was reported in Ref. 5, for 
both so*t inplane and stiff inplane rotors. 
This verifies the symbolic and numerical 
programs for the single-blade results and 
forms the basis for checking th <2 symboli- 
cally derived multiblade equations (and 
numerical results) subsequently. 

Figures 6 and 7 present the lead-lag 
damping plotted versus advance ratio from 
a flap-lag model, flap-lag-torsion model 
(two modes for each degree of freedom) , and 
flap-lag-torsion model with only one tor- 
sion mode, for a soft inplane rotor 
vWy * 0.7). Figure 6 presents the damping 
results for full elastic coupling (R - 1,0). 
It can be seen that the flap-lag model 
underpredicts the lead-lag damping. The 
model with only one torsion mode increases 
the damping above that of the model with 
two modes each. Figure 7 shows the lead- 
lag damping value plotted for zero elastic 
coupling (R * 0.0), The damping levels, are 
very m* ch reduced compared with those in 
the fuj.1 elastic coupling case. However, 
the flap-lag model is again the least 
damped. 

Lead-lag dampi g is plotted for a 
stiff inplane rotor (w v * 1.4) with varying 
advance ratio in Figs. 8 and 9. Figure 8 
presents the damping results for full 
elastic coupling. The same trend that was 
observed for the case of soft inplane 
(Fig. 6) exists. Here, it is to be noted 
that root splitting for high advance ratios 
occurs even when only ore torsion mode is 
used. Figure 9 presents the lead-lag damp- 
ing for increasing advance ratio for a 
stiff inplane rotor for zero elastic cou- 
pling parameter. Although a flap-lag model 
predicts a stable system, the rotor is 
unstable. This demonstrates the importance 
of elastic blade torsion in a forward- 
flight stability analysis. 

Figure 10 shows the lead-lag damping 
plotted for an advance ratio of y ■ 0.25 


while varying the elastic coupling param- 
eter for a stiff inplane rotor. Here a 
flap-lag model predicts positive damping 
for all values of R, whereas for a flap- 
lag-torsion model the damping varies with 
the elastic coupling value, increasing 
with the elastic coupling parameter. 

Multiblade-Equation Results 

The following figures present the 
lead-lag regressing mode damping results 
obtained from multiblade equations. The 
multiblade equations were explicitly 
derived using the symbolic program. This 
required explicit definition of all non- 
linear contributions and degrees of free- 
dom in terms of their harmonics. The 
result was a significant increase in the 
amount of data required by the symbolic 
program. Since this is a feasibility 
study on the use of symbolic programs in 
FORTRAN, only first harmonics were con- 
sidered in the nonlinear contributions. 
Consequently , damping data determined from 
the multiblade equations may differ from 
the single-blade solution. Additionally, 
the multiblade results are obtained by 
retaining only one torsion mode, although 
the nonlinear contribution from both tor- 
sion modes is used. This significantly 
reduces the time required ior the Floquet 
stability analysis. 

The damping values were first checked 
with those obtained from a single-blade 
solution obtained previously to validate 
the multiblade equation derivation pro- 
cess. It was found that the approximation 
cos B * 1.0, used in deriving the explicit 
multiblade equationc, will predict slightly 
higher (but less than 2%) damping for stiff 
inplane rotors with full elastic coupling 
parameter; this is because the approxima- 
tion has its greatest effect on the cou- 
pling elements. For all other values of 
the elastic coupling parameter, this 
approximation does not affect the resul- 
tant damping values. Where required for 
comparisons, the single-blade damping 
values are recalculated using this approx- 
imation; this is done to avoid the rederi- 
vation of the multiblade equations. 

Figures 11 and 12 show the lead-lag 
regressing mode damping plotted for a 
varying advance ratio with and without 
dynamic inflow from a flap-lag-torsion and 
flap-lag model for a soft inplane rotor. 
Figure 11 shows the damping for full elas- 
tic coupling (R ■ 1.0). For the flap-lag- 
torsion model, the dynamic inflow reduces 
the damping at practically all advance 
ratios. Its effect is negligible at 
advance ratios of 0.15 to 0.25. for the 
flap-lag model, the dynamic inflow 
increased the damping up to an advance 
ratio of 0,33; at higher advance ratios it 
reduced the damping. Figure 12 presents 
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the damping results for the zero elastic 
coupling (R = 0,0) . It is seen that for 
the flap-lag-torsion model, the dynamic 
inflow again reduced the damping in hover. 
Vet at intermediate advance ratios, dynamic 
inflow increased the damping, and at higher 
advance ratios, it once again reduced the 
damping. For the flap-lag model, the 
dynamic inflow increased lead-lag damping 
for all advance ratios. This is consis- 
tent with the flap- lag model results of 
previous studies (e.g.. Ref. 7). 

The lead- lag regressing mode damping 
is plotted for a stiff inplane rotor for a 
varying advance ratio in Figs. 13 and 14. 
Figure 13 is for a rotor with full elastic 
coupling. For a flap-lag-torsion model, 
the dynamic inflow reduced the damping up 
to an advance ratio p of 0.41. For 
p > 0.41, this model slowed a slightly 
increased damping value. The flap-lag 
model with dynamic inflow shows a small 
increase in deunping. This damping incre- 
ment gets smaller with hxgher advance 
ratios. Figure 14 is for a rotor with 
zero elastic coupling. For this config- 
uration, the dynamic inflow increases 
damping for all advance ratios. Conse- 
quently both the f lap-lag-torsion and 
flap-lag model show the same trend. 

Figure 15 shows the lead-lag regress- 
ing mode damping plotted for a stiff 
inplane rotor at an advance ratio of 0.25, 
for varying elastic coupling. For the 
flap-lag-torsion model, dynamic inflow 
reduces the damping for R > 0.3, but it 
increases the available damping for 
R < 0.3. However, this increase is not 
sufficient to stabilize the inplane mode. 
With the flap-lag model, dynamic inflow 
shows an increase in damping for all valuis 
of elastic coupling. This is the same 
trend as was observed in Ref. 11 for the. 
case of hover with both flap- lag and flap- 
lag-torsion blade models. 

Constant-Coefficient Approximation Results 

The effect of a constant-coefficient 
approximation (CCA) is presented in Fig. 16, 
whore the real part of the exponent is 
plotted for a stiff inplane rotor with full 
elastic coupling ttu.c did show a splitting 
of roots with a full periodic coefficient 
analysis (Fig. 13) . The CCA does not show 
this splitting, since the frequencies are 
very much away from the real axis. For 
this analysis, the regressing and collec- 
tive modes did predict the same damping 
trend with advance ratio as shown by the 
full Floquet analysis. However, the pro- 
gressing mode showed poor agreement between 
a CCA analysis and a Floquet-theory analy- 
sis. This is because the constant- 
coefficient approximation will only be 
good for low-frequency modes. 


Conclusions 

A symbolic manipulation program 
written in FORTRAN was used to derive the 
aeroelastic analysis equations of an 
elastic blade with flap-lag-torsion 
degrees of freedom in forward flight. 

The feasibility o- using the program to 
obtain explicit equations in a harmonic 
balance method and multiblade equations 
was studied. Numerical results were pre- 
sented, with and without dynamic inflow 
for a propulsive trimmed rotor. Both a 
flap- lag-tors ion model and a flap- lag 
model were analyzed. Soft inplane and 
stiff inplane rotors were considered. 

The following conclusions were drawn 
from this study of the use of a symbolic 
program for predicting rotor aeroelastic 
stability. 

1) The symbolic program can be used 
to obtain explicit equations. 

2) With the present program capa- 
bility, tne amount of data to the sym- 
bolic program increases greatly with the 
number of harmonics and degrees of 
freedom, 

3) In deriving the explicit harmonic 

balance equations and multiblade equations, 
the following should be noted: a) to 

obtain the harmonic balance equations, a 
numerical method is suggested since an 
arbitrary number of harmonics can be used 
without increasing the input data to the 
symbolic program; b) to obtain the multi- 
blade equations, the perturbed equations 

in their Fourier series form are derived 
using the symbolic program. Then the 
multiblade equations themselves are 
obtained numerically. 

It is recommended that a selective 
(judicious) combination of symbolic and 
numerical programs is required for an 
efficient derivation and numerical -study 
process. 

The following conclusions were drawn 
from the numerical study of a single-blade 
solution. 

1) A flap- lag-torsion stability 
analysis from a trim procedure in which 
only the flap degree of freedom is used 
underpredicts the lead-lag damping. 

2) In the case of stiff inplane 
rotors, high forward flight speed is 
destabilizing. At high advance ratios, a 
splitting of the roots is encountered, 
yielding two real-part characteristic 
exponents a' '.he same frequency. 
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3) Using only one torsion mod* usu- 
ally increases the damping value f.om the 
flap-lag structural model. 

4) The damping values for a stiff 
inplane rotor are very sensitive to elas- 
tic coupling parameter Depending on this 
parameter , the rotor ctn be either stable 
or unstable. 

The following conclusions were drawn 
from the numerical study of a multiblade 
solution with dynamic inflow. 

1) For a flap-lag model, and for both 
soft inplane and stiff inplane rotors with 
zero elastic coupling, the dynamic inflow 
increased damping at all advance ratios 
considered; with full elastic coupling, the 
dynamic inflow increased the damping at low 
advance ratios, but reduced damping at high 
advance ratios. 

2) For a flap-lag-torsion model, 
dynamic inflow slightly reduced lead- lag 
regress ing-mode damping for full elastic 
coupling. The same trend was observed for 
both soft inplane and stiff inplane rotors. 

3) For a given advance ratio, the 
variation of damping with elastic coupling 
parameter for a stiff inplane rotor showed 
the same trend as did the hover case. 

4) The constant-coefficient approxi- 
mation for the stiff inplane rotor does not 
show the splitting of the roots, since the 
frequency of the lag mode is away from the 
real axis. 
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Table 1. FORTRAN symbols 


Original 

variable 

FORTRAN 

symbol 

Original 

variable 

FORTRAN 

symbol 

R 

RAD 

v' 

VS 

u 

U 

Vi 

LAMB 

u 

UD 

w 

W 

Ut 

UT 

w' 

ws 

Up 

UP 

w 

WD 

Vf 

UF 

X 

XCOR 

Up 

UFD 

0 

THTA 

V 

V 

4> 

PHI 

V 

VD 

n 

OMEG 



Table 2. Typical input to HESL and output to 
calculate tangential and normal velocities 
{ U«ji and Up) 


READ MATRIX 

ILAFP0303 

03 

1.0 

-0.5 VS VS 
-0.5 WS WS 
01 

1.0 VS 

01 

1.0 WS 

02 

-1.0 VSIEONE 

-1.0 WStETWO 

03 

1 . 0 tETRE 
-1.0 PHI THTA 
-1.0 VS WS THTA 
02 

1 . 0 PHI 
1 . 0 %EFUR THTA 
02 

1.0 VSSETWO 

-1.0 WStEONE 

03 

-1.0 %ETRE THTA 
-1.0 PHI 

-1.0 VS WS 


- 02 

-1.0 PHI THTA 

1 . 0 %EFUR 
READ MATRIX 
IVEL0301 
04 

-1.0 UD 

1.0 UFD 

1.0 OMEG V 

1.0 MU OMEG RAD CSCY 

03 

-1.0 VD 

-1.0 tXUUF OMEG 

-1.0 MU OMEG RAD SNCY 

02 

-1.0 OMEG RADtLAMD 

-1.0 WD 

FORM MATRIX 

ILAFF t VEL#AVEL*E2D1 

MATRIX EXPRESSION 

02IAVEL 

% UT0201* UP0301 

NEGATE EXPRESSION 
% UT 

NEGATE EXPRESSION 
% UP 


Note: %EONE, %ETWO etc. are expressions read 

earlier in the program. 


Output of Up 


•A*******.********* ******** **************** 

* DETAILS OF THE EXPRESSION UT * 

* NUMBER OF TERMS 22 * 

******************************************* 


1 

1.000* 

VS*OMEG* V* 


2 

1.000* 

VS*0MEG* MU* RAD*CSCY* 

3 

1 . 000*THTA* WS*0MEG* MU* 

RAD*CSCY* 

4 

1.000* 

PHI* WS*0MEG* MU* 

RAD*CSCY* 

5 

1.000* 

VD* 


6 

1 . 000*XCOR*OMEG* 


7 

1.000* 

U*OMEG* 


8 

-1.000* 

UF*OMEG* 


9 

1.000*OMEG* MU* RAD*SNCY* 


10 

-0.500* 

VS* VS* VD* 


11 

-0.500* 

VS* VS*XCOR*OMEG* 


12 

-0.500* 

VS* VS*OMEG* MU* 

RAD*SNCY* 

13 

-0,500* 

PHI* PHI* VD* 


14 

-0.500* 

PHI* PHI*XCOR*OMEG* 


15 

-0.500* 

PHI* PHI*OMEG* MU* 

RAD*SNCY* 

16 

-1,000* 

PHI*THTA* VD* 


17 

-1,000* 

PHI*THTA*XCOR*OMEG* 


18 

-1.00Q* 

PHI *THTA*0MEG* MU* 

RAD*SNCY* 

19 

1.000* 

PHI*LAMB*OMEG* RAD* 


20 

1,000* 

PHI* WD* 


21 
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AXIS OF ROTATION 


0n 



Fig. 1 Rotor-blade coordinate systems and 
deflections. 



Fig. 3 Flowchart of the aeroelastic 
analysis. 


; 



Fig. 2 Inflow dynamics. 
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Fig. 4 The effect of the number of 
degrees of freedom used in trim analysis 
on lead-lag damping versus advance ratio: 
soft inplane, u» v * 0.7, R = 1.0, propulsive 
trim. 
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FLAP-LAC-TORSION TRIM 
FLAP TRIM 



NUMBER OF MODES 
o 2 FLAP, 2 LAG, 2 TORSION 
a 2 FLAP, 2 LAG, 1 TORSION 
□ 2 FLAP, 2 LAG 


0 .05 .10 .15 .20 ^25 !3o ^35 40* 45 

ADVANCE RATIO, u 

Fig. 5 The effect of number of degrees of 
freedom used in trim analysis on lead-iag 
damping versus advance ratio: stiff 9 

trim?" 6 * ” 1 ’ 4 ' *-1.0. propulsive 



•06 .10 .15 .20 .25 .30 .35 40 

“■ ADVANCE RATIO, u 

Fig. 7 Lead-lag mode damping versus 
«d a ?lf^° f ° r , a f l»P-lag-torsion model 

r.V 8P ; la V° dels soft in Plane f 

“V 0.7, R * 0.0. 


-.06 


K -.016 


NUMBER OF MOOES 
o 2 FLAP ■> LAG, 2 TORSION 
c 2 FLAP, 2 LAG, 1 TORSION 



■ 10 .15 .20 .25 .30 

ADVANCE RATIO, ju 

Fig. 6 Lead-lag mode damping versus 

and** 0 ? 1 rat f° for a f lap-lag-torsion model 
and a flap-lag model: soft inplane, 

« v " 0.7, R - l.o. 


•06 .10 .16 .20 .26 .30 36 
ADVANCE RATIO, *i 

Fig. 8 Lead-lag mode damping versus 

and** 0 ?! ^1° for . a . £ l*P-lag-torsion model 
and a flap-lag model: stiff inplane, 

Wy - 1.4, R - l.o. 
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NUMBER OF MODES 
G 2 FLAP. 2 LAG. 2 TORSION 
" L 2 FLAP. 2 LAG. 1 TORSION 
□ 2 FLAP. 2 LAG 

2ND LEAD-LAG MODE. FLAP-LAG-TORSION 

(2 MODES EACH) 


NUMBER OF MODES 
o 2 FLAP. 2 LAG. 2 TORSION 
a 2 FLAP. 2 LAG. 1 TORSION 



Fig. 9 Lead- lag mode damping versus 
advance ratio for a flap-lag-torsion model 
and a flap-lag model: stiff inplane, 

cjy * 1.4, R = 0.0. 
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Fig. 10 Lead-lag mode damping versus 
elastic coupling for a flap-lag-torsion 
model and a flap-lag model: stiff inplane, 

u) v = 1.4, y = 0.25. 







Fig. 11 The effect of torsion and dynamic 
inflow on lead-lag regressing mode damping 
versus advance ratio: soft inplane, 

cj v = 0.7, R - 1.0 (multiblade equations). 


o FLAP-LAG TORSION 



Fig. 13 The effect of torsion and dynamic 
inflow on lead-lag regressing mode damping 
versus advance ratio: stiff inplane, 

w v = 1.4, R * 1.0 (multiblade equations). 
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Fig, 12 The effect of torsion and dynamic 
inflow on lead-lag regressing mode damping 
versus advance ratio: soft inplane, 

(o v «0,7, R ■ 0.0 (multiblade equations). 
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Fig. 14 The effect of torsion and dynamic 
inflow on lead-lag regressing mode damping 
versus advance ratio: stiff inplane, 

u) v -1.4, R - 0,0 (multiblade equations). 
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Fig. 15 The effect of torsion and dynamic Fig. 16 Comparison of constant-coefficient 
inflow on lead-lag regressing mode damping approximation and Floquet analysis: stiff 

versus elastic coupling, stiff inplane, inplane, w v = 1.4, R = 1.0 (multiblade 

u> v ■ 1.4, \i = 0.25 (multiblade equations). equations). 


i 



i 




DISCUSSION 
Paper Ih. 15 



THE INFLUENCE OF DYNAMIC INFLOW AND TORSIONAL FLEXIBILITY ON ROTOR DAMPING IN FORWARD FLIGHT 

FROM SYMBOLICALLY GENERATED EQUATIONS 
T, S. R. Reddy 
and 

William Warmbrodt 


Peretz Friedmann, University of California, Los Angeles : I'd like to congratulate you on a very 

nice paper. Obviously I have a vested interest because for the last three years I have been 
waiting for somebody to redo the problem to find out whether Kottapalli and I have done it 
correctly. Now that you have shown these results and Neelakanthan has shown some results at the 
last European Forum where also the same trends were exhibited I guess I can sleep in peace. 

What I really wanted to emphasize are two things. One is the contribution you have made is a 
really significant one because as somebody who has derived equat »*s by hand for a long time I 
definitely believe that the way to go is to use a computer. The ond comment which I have— 
and it is in the form of a question— is if I correctly understand cne results you have shown 
then it seems to be that dynamic inflow doesn't have an awful lot of influence in the case of 
coupled flap-lag-torsion in forward flight. I was wondering if you would agree with this 
statement? 

Reddy : Yes, that's what our results show. 

Wayne Johnson. NASA Ames Research Center : Following along these lines of what Peretz was dis- 
cussing about using the computer: in dealing with this subject, if you were going to do the 

same work, but do it over again were there any pieces of the problem that you did with pencil 
and r^per, [that you] did by hand that if you had to do it all over again you would automate 
these pieces also. In other words, was there anything left to put into the computer? 

Reddy : Yes, we are now finding that the required input data has increased tremendously so we 

will have to change the program to minimize data Inputs. 

Friedrich Straub. Hughes Helicopter s: How long did it take you to include the forward flight in 

the equations coming from the hover results? 

Reddy : I attended the SDM Conference on the 14th to the 16th of May. Then we came back to this 

area and submitted the abstract— that was the end of May. We finished the paper by August, but 
most of this time was spent on developing uhe trim and response solution program, so it took 
even less time. 

Bill Warmbrodt. NASA Ames Research Center : You might point out that the original derivation was 
done including the influence of forward flight; uowever, the program was first exercised to 
develop the hover results presented at the SDM Conference. 
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